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Abstract 

CFTR is a dynamically regulated anion channel. Intracellular WNKl -SPAK activation causes CFTR to change pernneability and 
conductance characteristics from a chloride-preferring to bicarbonate-preferring channel through unknown mechanisms. 
Two severe CFTR mutations (CFTR^'^^) cause complete loss of CFTR function and result in cystic fibrosis (CF), a severe genetic 
disorder affecting sweat glands, nasal sinuses, lungs, pancreas, liver, intestines, and male reproductive system. We 
hypothesize that those CFTR mutations that disrupt the WNKl -SPAK activation mechanisms cause a selective, bicarbonate 
defect in channel function (CFTR^'^) affecting organs that utilize CFTR for bicarbonate secretion (e.g. the pancreas, nasal 
sinus, vas deferens) but do not cause typical CF. To understand the structural and functional requirements of the CFTR 
bicarbonate-preferring channel, we (a) screened 984 well-phenotyped pancreatitis cases for candidate CFTR*" mutations 
from among 81 previously described CFTR variants; (b) conducted electrophysiology studies on clones of variants found in 
pancreatitis but not CF; (c) computationally constructed a new, complete structural model of CFTR for molecular dynamics 
simulation of wild-type and mutant variants; and (d) tested the newly defined Cfr/?*° variants for disease in non-pancreas 
organs utilizing CFTR for bicarbonate secretion. Nine variants {CFTR R74Q, R75Q, Rn7H, R170H, L967S, L997F, D1152H, 
S1235R, and D1270N) not associated with typical CF were associated with pancreatitis (OR 1.5, p = 0.002). Clones expressed 
in HEK 293T cells had normal chloride but not bicarbonate permeability and conductance with WNKl -SPAK activation. 
Molecular dynamics simulations suggest physical restriction of the CFTR channel and altered dynamic channel regulation. 
Comparing pancreatitis patients and controls, CFTR*" increased risk for rhinosinusitis (OR 2.3, p<0.005) and male infertility 
(OR 395, p<<0.0001). WNKl -SPAK pathway-activated increases in CFTR bicarbonate permeability are altered by CFTR^° 
variants through multiple mechanisms. CFTR^'^ variants are associated with clinically significant disorders of the pancreas, 
sinuses, and male reproductive system. 
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Introduction 

The cystic fibrosis transmembrane conductance regulator 
(CFTR, GenBank Accession: AH006034.1) is an ATP-binding 
cassette (ABC) transporter-type protein localized to the apical 
plasma membrane of epitheKal cells. It difTers from other ABC 
transporters in that it acts as a regulated anion channel rather than 
a transporter [1]. When the channel is open, anions move across 
the membrane down their electrochemical potential gradient, 
resulting in fluid and electrolyte secretion or absorption. 

The CFTR molecule has been intensely studied because 
mutations in the CFTR gene are associated with cystic fibrosis 
(CF, OMIM #219700), the most common life-threatening genetic 
disorder among populations of Northern European ancestry [2, 3]. 
However, the clinical features of CF and CFTR-related disorders 
are variable, and laboratory studies of CFTR regulation, its 
biophysical properties and molecular mechanisms of (dys)function 
have been challenging due to the complexity of the regulatory 
mechanisms and the dynamic flexibility of various structural 
domains (see recent reviews [4,5]). 

Cystic fibrosis is an autosomal recessive syndrome usually 
caused by inheriting two CFTR mutations that eliminate effective 
chloride conductance {CFTR^'')[2,3]. Although nearly 2000 
CFTR variants have been described (http://www.genet.sickkids. 
on.ca), the majority of CF cases are associated the CFTR 508F-del 
mutation as a homozygous genotype or in combination with 
another severe CF-associated mutation {CFTR'^''' /CFTR'^''') that 
together result in minimal CFTR function. Thus, most research 
has focused on the regulation of chloride conductance, and 
dynamic modeling of the first of two nucleotide-binding domains 
(NBDl), which normally contains F508 [4,5]. Based on numerous 
studies, three conformations have been described for the molecule 
as an anion channel: a closed state, an open state, and an open- 
ready state [4]. However, the relative permeability/conductance 
ratios of chloride and bicarbonate are variable [6] and may be 
dynamically regulated [7] , suggesting that conformational changes 
induced by point mutations in the channel or in the permeability 
pore may alter ion permeation properties of CFTR. 

Diagnosis of CF is based on a combination of phenotypic 
features, family history, functional tests and/ or the identification 
CFTR^'^ variants on both alleles [8,9]. Organ dysfunctions begin 
in utero and include chronic pancreatitis, meconium ileus, and 
congenital bilateral agenesis of the vas deferens. Progressive 
sinorespiratory dysfunction develops in childhood due to bacterial 
infections, inflammation, and scarring, and male infertility is 
recognized in adulthood. Disease severity and complexity is 
modified by other genes [10-1 2], environmental factors [13], and 
mild-variable CFTR variants [3,14]. Mild CF phenotypes, CFTR- 
related disorders limited to a single organ, are associated with non- 
CFTR^^ variants with residual channel function, classified as 
mHd-variable variants {CFTR"'-) [2,9,15]. 

CFTR'^'^ and CFTR'"'^' variants are associated with recurrent 
acute pancreatitis and chronic pancreatitis [16-19]. Recently, we 
reported that the variant CFTR R75Q, which was previously 
classified as benign, is associated with familial and sporadic 
chronic pancreatitis, either with another CFTR variant (recessive) 
or with the serine protease inhibitor, Kazal Type 1 {SPINKl) 



N34S high-risk haplotype (complex genotype) [18]. Patch-clamp 
studies of CFTR R75Q^ clones under standard conditions 
demonstrated normal chloride conductance but a selective 
disruption in bicarbonate conductance [18]. Thus, CFTR R75Q 
causes selective bicarbonate defective {CFTR^^) conductance and 
is associated with chronic pancreatitis but not CF[18]. It is not 
known if other CFTR variants share this phenotypic feature, 
whether the defect is associated with the channel function under 
all or special conditions, or if other mechanism(s) underlying these 
observation. 

Independendy, we demonstrated that CFTR bicarbonate 
(HCO3 ) permeability increases through WNKl-SPAK signaling 
pathway activation [20]. WNKl is member of the "with-no-K" 
(Lys) kinases that serves as a sensor of osmolality, chloride 
concentration, and other factors within cells and respond by 
activating additional kinases linked to a variety of ion channels and 
exchanges, including CFTR [20-22]. In cell-based models, low 
intracellular chloride concentrations ([CI"];) result in WNKl- 
mediated SPAK activation that strongly increases CFTR HCO3 
permeability in CFTR-transfected HEK 293T, PANCl, and 
guinea pig pancreatic duct cells, making CFTR primarily an 
HCO3 channel [20]. The structural and dynamic mechanisms of 
this phenomenon are unknown. 

We hypothesized that CFTR variants that disrupt the WNKl- 
SPAK-associated increase in bicarbonate permeability wiU in- 
crease the risk of pancreatitis and affect other organs in which 
CFTR is used for bicarbonate secretion. To test this hypothesis 
and to gain insight into potential mechanisms, we adopted a 
multidisciplinary approach. First, to identify candidate CFTR^^ 
variants, we conducted a systematic review of the literature to 
compile CFTR variants that have been reported at least twice in 
previous chronic pancreatitis case-control genetic studies, plus 
common CFTR'"'' variants. Second, using this panel of 81 CFTR 
variants (Table SI in the Supplementary Material), we genotyped 
the deeply phenotyped North American Pancreatitis Study 2 
(NAPS2) subjects[23] to identify candidate CFTR'^° variants that 
were also present in our cases and controls (43 of them, listed in 
Table 1). Third, to determine if CFTR'^^ variants are associated 
with altered WNKl-SPAK pathway-stimulated CFTR bicarbon- 
ate permeability, we generated plasmids containing the candidate 
CFTR'^^ variants selected from the NAPS2 study and expressed 
them in HEK-293T cells for electrophysiological analysis. Fourth, 
to gain insight into the molecular mechanisms of CFTR^^ 
dysfunction, we performed molecular dynamics (MD) simulations 
based on homology-modeled structures of ABC transporters and 
examined the effect of CFTR^^ variants on the structure and 
dynamics of the channel. Finally, to determine if CFTR^^ variants 
are associated with disease in non-pancreatic tissues, we used the 
phenotyping criteria for sinusitis and male infertility for the 
NAPS2 cases and controls. 

These studies revealed at least 9 CFTR variants. We found 
that the WNKl-SPAK pathway that enhances CFTR bicarbonate 
permeability/conductance compared with chloride conductance 
in HEK-293T ceUs is altered by CFTR'^^ variants. The 
examination of MD trajectories suggests at least two potential 
mechanisms of channel dysfunction. Phenotype-genotype studies 



PLOS Genetics | www.plosgenetics.org 



2 



July 2014 I Volume 10 | Issue 7 | e1004376 



CFTR Bicarbonate Permeation Variants 



Author Summary 

Genetic disorders of ion channels can affect the body's 
ability to function properly in many ways. CFTR, an ion 
channel regulating movement of chloride and bicarbonate 
across cell membranes, is important for absorbing and 
secreting fluids. If the gene responsible for the CFTR 
channel is mutated severely, the result is cystic fibrosis, a 
hereditary disorder in which the patient develops thick 
mucus, especially in the lungs, as well as scarring (fibrosis) 
in the pancreas. Cystic fibrosis also affects the sweat 
glands, nasal sinuses, intestines, liver, and male reproduc- 
tive system. Mutations to the CFTR gene that do not cause 
cystic fibrosis have been considered benign. However, we 
discovered 9 CFTR mutations that do not cause cystic 
fibrosis but do cause inflammation and scarring of the 
pancreas (chronic pancreatitis). These mutant CFTR chan- 
nels secrete chloride, which is important in the sweat 
glands, lungs, and intestines, but not bicarbonate, which is 
important in the pancreas, sinuses, and male reproductive 
tract. We found patients with any of these 9 mutations had 
chronic pancreatitis, and often sinus infections, and male 
infertility, but not other symptoms of cystic fibrosis. Our 
computer models and data will help researchers develop 
better drugs and help physicians treating patients with 
chronic pancreatitis. 

in humans demonstrated that CFTR^^ variants are also associated 
with disorders of the pancreas, sinuses, and male reproductive 
systems. 

Results 

CFTR genotyping in pancreatitis patients and controls 

We genotyped 984 well-phenotyped cases of pancreatitis from 
NAPS2 for 81 CFTR variants, including common CF mutations 
and variants previously reported in at least two subjects with 
pancreatitis but not CF. Common tag-SNPs at the CFTR locus 
were previously excluded in a pancreatitis genome-wide associa- 
tion study (all p values ^0.01) [24], suggesting that the missing 
heritabUity and predicted dysfunction was primarily associated 
with multiple rare variants. SPINKl N34S was also genotyped to 
determine complex risk [18]. Only SPINKl N34S heterozygotes 
were used for trans-heterozygote analysis with CFTR, since 
homozygous SPINKl N34S is sufficient to cause pancreatitis. 

Of 43 CFTR variants identified in the NAPS2 cohort (Table 1), 
nine not associated with typical CF but reported in patients with 
pancreatitis[25-29] were of particular interest: R74Q, R75Q, 
R117H {CFTR'"-'' only when m cis with IVS8-T5[30]; 
R117H*T5), R170H, L967S, L997F, D1152H, S1235R, and 
D1270N. These were either independendy associated with disease, 
were found in subjects with SPINKl N34S as a complex high-risk 
trans-heterozygous genotype or had predicted clinical relevance 
based on prior reports or their location on the CFTR molecule. 
Taken together, these nine CFTR^^ variants were found more 
commonly in cases (14.2%) than controls (9.8%) (OR 1.5, 
p = 0.002) (Table 1). 

As expected, CFTR variants associated with typical CF were 
also identified in more cases than controls (8.7% cases, 3.3% 
controls; OR 2.8, p<0.0001). Other candidate CFTR variants, 
including I148T, M470V, T854T, Cll463Q,and the "5T" aUele, 
were either rare or were not associated with pancreatitis in our 
cohort (Table 1). A total of 189 cases (19.8%) carried one or more 
CFTR variants of any kind (controls 13.0%, p<0.0001, OR 1.6, 



95% C.I 1.3-2.0): 38% of these mutations were CFTi?'^'''^ variants, 
while the remaining were 

QpjjlBD {&2%). Twenty-five 

cases and no controls carried multiple mutations in CFTR. 
Twenty-five cases carried trans-heterozygous mutations in both 
CFTR and SPINKl (N34S), including five patients with three or 
more mutations (Table 2). 

Several candidates that were previously reported to be 
associated with pancreatitis or atypical CF were not replicated in 
the NAPS2 cohort. I148T was seen in three cases and one control, 
so an efTect could not be detected or excluded; the in cis deletion 
mutation 3199del6 was not detected in any I148T carriers. The 
IVS8T5 variant was identified in 9.9% of cases and 8.2% of 
controls, which is not individually significant. There were six 
N34S/T5 trans-heterozygote controls and no cases, but the 
combined efTect of the SPINKl N34S variant with IVS8T5 was 
not significantly higher than N34S alone. Four variants were 
identified in only one patient and no controls: CF mutations 
2184delA, 3120+lOA, R1162X, and mutation of varying 
clinical consequence, G1069R. 

Functional assays on CFTR variants 

For our functional studies, we cloned the nine CFTR variants 
and confirmed that they had normal folding, glycosylation 
(Figure la) and chloride channel activities, except for R117H 
(Figure lb). Because CFTR bicarbonate permeability is dynam- 
ically increased through [CI ]i-sensitive WNKl-SPAK signaling 
pathway activation [20], we tested this in HEK 293T cells[20] 
using whole-cell current measurements by replacing 150 mM 
extracellular CF with 140 mM HCO3" and 10 mM CF. 
Representative traces for voltage and current measurements are 
presented in Figures Ic, Id, and S2, and a summary of the 
indicated numbers of experiments is depicted in Figure le and f. 

The bicarbonate permeability of CFTR in cells that do not 
overexpress WNKl and SPAK was much smaller than that of 
chloride, with Phcos/Pci = 0.24±0.05 (Figure SI). As reported 
previously [20], with WNKl and SPAK co-expression and low 
[CF];, the permeability of CFTR to bicarbonate increased and 
reached that of chloride, with Phco:vPci= 1.06±0.06 (Figure Ic 
and le). In contrast, CFTR Phco3/Pgi failed to increase in 
CFTR R170H (Figure Id) and all of die candidate CFTR'^'^ 
variants (Figures le and S2). Furthermore, all CFTR^^ 
candidate variants lowered bicarbonate conductance (Ghco:}/ 
Oca), which is an important metric determining apical bicarbonate 
efflux in CFTR-expressing epithelia (Figure If); the decrease was 
statistically significant for all variants except D1270N. Treatment 
with the CFTR inhibitor CFTR;„h-172 (20 |xM) inhibited >90% 
of the HCOs" currents (Figure S2), indicating that CFTR 
mediates most of the HCO3 currents observed in the present 
experiments. 

To further evaluate the mechanism of bicarbonate conductance, 
we tested the hypothesis that the well-established CFTR channel 
blocker, CFTR;„h-172[31-33] blocks HCO3" current. We found 
that CFTR;„h-172 (20 |J,M) inhibited >90% of the HCO3" 
currents (Figure S2), indicating that with WNKl-SPAK activation, 
CFTR mediates most of the HC03~ currents. 

Structural and dynamic modeling of CFTR wild-type and 
variants 

The specific amino acid substitutions that interfere with 
WNKl/SPAK-activated transformation of CFTR to a more 
efficient bicarbonate-conducting channel are scattered throughout 
the linear DNA sequence, suggesting that three-dimensional 
structure and/or mechanisms of dynamic conformational changes 
linked to these amino acids are important risk for pancreatitis. We 
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> computationally modeled the molecular structure, and studied 

^ o the dynamics, of wild type (WT) and mutated CFTR channels. 

^ S Because the efTective van der Waals radius of chloride (1.8 A 

1 £ [34]) is smaller than that of bicarbonate (2.6 A, see Methods), we 
° ^ tested whether amino acid substitutions that reduced the inner 
^ diameter of the CFTR channel could selectively impede 

2 J bicarbonate conductance. A CFTR-WT model (Figure 2a) 
S g was constructed [35,36] and used to locate and study CFTR'^^ 
5 S functional variants (Figure 1). The model is based on the most 

recendy resolved ABC transporter structure (from Staphylococcus 
aureus savl866; see Materials and Methods). Figure S3 a shows 
> o the superposition of our model on this crystal structure, which 

g u yields an RJvISD of 1.6 A. Panel b shows that residues lining the 

>- ^ pore at the membrane-spanning domain (MSD), observed by the 

5 3 end of 50 ns simulations, agree in general with the CFTR model 

^ J built by Norimatsu and collaborators [37,38] which was also 

£ "g confirmed by cysteine scanning experiments [37,39]. Likewise, 

g o the pore radius profile evaluated for our wild-type structural 

§ ^ model (Figure S3 c solid curve, with the gray band displaying 

™ the fluctuations observed in 50 ns simulations) is qualitatively 

I 00 consistent with that observed by Norimatsu and coworkers [37] 

£ f for the MSD. 

3 s MD simulations comparing the channel diameters of the WT 
■-. X and mutants L997F and D1152H (Figure 2c— f) demonstrate 
^ 5" that the channel diameter is observed to narrow down from an 
c "g average value of 10.3 A to 7.5 A (standard deviation, ct = 0.5 A) 

at the pore region, near the L997F amino acid substitution 
.1 "g (Figure 2e), and from an average of 9.9 A to 4.3 A (a = 1. 1 A) 

g 1^ for the CFTR'^^ mutant D1152H (Figure 2f). Note that in 

contrast to the WT CFTR and L997F mutant where the 
structure maintains its stability, the D1152H mutation induces 
significant fluctuations in local conformation, which are reflected 
on the changes in the pore diameter at this location within the 
^ 5 ^ channel. 

Q i§ > In order to determine residues that play a key role in the global 

"5 § §■ dynamics of the CFTR, we performed an elastic network model 

c j= S (ENM) analysis. ENM analysis provides information on the 

" S z mechanisms of collective movements intrinsically accessible to the 

^ c z structure, which usually enable structural changes relevant to 

S * function [40]. Application to CFTR highlighted the critical 

positioning of R74, R75, R170, L967, and R1162 at the hinge 
region that modulates the collective movements of the nucleotide- 
binding domains (NBDs) with respect to membrane-spanning 
domains (MSDs) (mode 1 in Figure 3). We also note that 
L967, L997, Dl 152, and Rl 162 act as anchors in collective mode 
2. In this mode, the two NBDs are observed to move in opposite 
directions (see color-code diagram in Figure 3). The relative 
movements of the two NBDs, is known to control channel gating, 
hence the significance of this mode, or the alterations in mode 2 
^ S § potentially caused by substitutions at the corresponding hinge 

o _5 _c site. 

S _^ "S These two results suggest that substitutions of amino acids (or 

" o 5 their side chains) at those particular regions could have an impact 
(u g 



fc "5 K 

£ Is 5 



iS + * 

^ o = 
+^ o ^ 



5 .£ "S 



on the collective dynamics of CFTR, and interfere with concerted 



oj o S m movements that would otherwise facilitate anion permeation. We 

" j= o 2 noted that the mean-square fluctuations in our model are minimal 

I S S £ at those particular residues (Figure S4), suggesting that mutations 

S I £ at those sites could not be accommodated without affecting the 

E s E E overall transporter structure and dynamics. Minimal mobility at 
S g- t= .2 
O — ^ 

c 5 o the NBD2 (Figure 3), and had weaker functional effects than 

5 = 1- T3 " 



those mutation sites originates from the contribution of global 



(most collective) modes. In contrast, the CFTR'^^ candidate 
Si o 5 variants D1270 and SI 235 are in close proximity on the surface of 



other CFTR variants (Figure 1). 
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Table 2. CFTR variants in subjects with chronic rhinosinusitis or male infertility (age >30 years). 





Rhinoslnusitis 


Yes 


No 


p-value 


OR 


CI 


Controls 


53 (10.2%) 


468 








Cases (all) 


151 (15.9%) 


798 


0.002 


1.67 


1.19-2.38 


0 CFTR"" 


111 


649 


0.021 


1.51 


1.05-2.18 


0 cnfP 


111 


649 


0.021 


1.51 


1.05-2.18 


1 cnfP 


14 


50 


0.011 


2.47 


1.18-4.91 


1 CFTR"" 


23 


78 


0.001 


2.60 


1 .43-4.60 


1 CFTR" or 1 CFTR"" 


37 


128 


0.0001 


2.55 


1.55-4.15 


CFTR'^/CFTR'" or CFTR''°/CFTR'^^ 


3 


21 


0.73 


1.26 


NS 


Infertility 


Yes 


No 


p-value 


OR 


CI 


Controls 


1 (0.4%) 


160 


_ 


_ 


_ 


Cases (all) 


1 7 (4.2%) 


390 


0.03 


7.14 


1.10-300 


0 CFTR"" 


8 


329 


0.28 


3.88 


0.5-174 


0 CFTR^'' 


8 


329 


0.28 


3.88 


0.5-174 


1 CFTlf 


2 


24 


0.051 


13.0 


0.65-786 


1 CFTR™ 


2 


35 


0.090 


8.99 


0.46-541 


1 CFIR^^ or 1 CFTR"" 


4 


59 


0.023 


10.7 


1.03-536 


CFTR'^/CFTR"" or CFTR^^/CFTR" 


5 


2 


1 .20E-07 


303 


23-15783 



Top: Chronic rhinosinusitis in NAPS2 controls and cases with 0, 1, or 2 CFTR mutations. 

Bottom. Self-reported prevalence of male infertility among males over 30 years of age. Odds ratios were calculated comparing CFTR carrier cases in each subcategory 
against all controls. Because CFTR'^^ and CFTR^" both affect bicarbonate conductance, we calculated the association and risk associated with the presence of either 
variant type (shaded). 
doi:l 0.1 371 /journal.pgen.l 004376.1002 



Association of CFTR variants with sinus disorders and 
male infertility 

To examine the potential clinical relevance of CFTR'^^ 
variants, we reviewed case report forms for additional CF 
phenotypic features of dysfunction in the sinorespiratory and male 
reproductive systems, which both use CFTR for bicarbonate 
secretion. Association with CFTR''''' alleles was used to test for 
CFTR-mediated chloride secretion, CFTR^'^ to test for selective 
bicarbonate-mediated secretion and, because both CFTR^'^ and 
CFTR'^'' cause defective bicarbonate conductance, association 
with either CFTR^'^ or CRTR^'^ alleles, or recessive genotypes 
{CFTR'^'^/CRTR'^'^ or CFTPS'''/CRTPv'^'^) to assess overall risk 
of altered bicarbonate secretion on organ dysfunction. 

The sinuses may use CFTR bicarbonate secretion, in part, for 
mucus hydration [41]. Sinusitis is common, with a complex gene- 
environment-anatomic risk that includes anatomy, allergies and 
recurrent infections. Self-reported chronic sinusitis was more 
common in pancreatitis cases (n= 151; 15.9%) than in controls 
(n = 53; 10.2%, P = 0.002) (Table 2). We identified die R75Q, 
Rl 17H, L967S, L997F, Dl 152H, and S1235R CFTR'^'' variants 
as well as Ci*'77?''^-associated variants (e.g., F508del, G542X) in 
cases with rhinosinusitis. Sinusitis was reported in pancreatic 
cases who did not have any of the CFTR variants in our test 
panel (p = 0.021; OR 1.51; CI 1.05-2.18), but risk increased 
among carriers of CFTi?^^ (p = 0.001; OR 2.60, CI 1.43-4.60), 
QfYjlCF ^ ^ Q Qj. 2.47; CI 1.18^.91) or either CFTR'^'^ or 
QpjllCF ^a^j.ia^nt allele (p = 0.0001; OR 2.55; CI 1.55-4.15) 
(Table 2). Rhinosinusitis was not statistically associated with 
recessive genotypes, possibly due to the complex nature of 
chronic sinusitis or requirement for an unidentified epistatic risk 
factor. 



CFTR bicarbonate secretion also plays a role in pH regulation 
in the male reproductive system [42]. Male infertility is uncommon 
and not dependent on recurrent infections. Self-reported male 
infertility over age 30 years was more common among cases 
(n= 17; 4.2%) than controls (n= 1; 0.4%, p = 0.03) (Table 2). We 
identified R75Q, R117H, and S1235R as weU as the CFTR'^'' 
variants F508del, G542X and 2789-H5G<A in male cases with 
infertility. There was no increased risk of male infertility in cases 
without CFTR variants (p = 0.28), but there was significant risk in 
cases with either CFTR'^° or CFTPS'' alleles (p = 0.023; OR 10.7; 
CI 1.03-536) or as a recessive genotype (p = 1.2 x 10"^; OR 303; 
CI 23-15783) (Table 2). 

Discussion 

Our integrative approach revealed a new functional class of rare 
CFTR variants of clinical significance in pancreatic disease. Targeted 
genotyping of reported and plausible CFTR variants in our cohort 
identified candidate variants with a high pre-test probability of being 
diseases associated, and these were evaluated for specific flinctional 
studies in model cell types and focusing on a context-dependent 
signaling pathway. Although the CFTi?^^ variants were scattered 
throughout the genetic sequence, three-dimensional models of the 
protein provided insight into structural and dynamic mechanisms of 
dysfiinction. Significant association between CFTR'^'^ variants and 
symptoms of sinusitis and male infertility, but not overt lung disease as 
in CF, provided additional evidence of context-depended dysfunction 
in humans. We believe that this type of integrated approach will be 
important in understanding the genetic contribution to this and other 
complex disorders and informing the development of therapeutics 
that target the molecular etiology rather than the phenotype. 
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Figure 1. Functional characteristics of the nine Cfr/f^" variants. Panel a. Wild-type (WT) and variant CFTR proteins were expressed in HEK 
293T cells and immunoblotted with anti-CFTR and anti-Aldolase antibodies. Replicate lanes are in small font. Band B, expected size of immature ER 
core-glycosylated CFTR; band C, mature complex-glycosylated CFTR. Panel b. Whole-cell Cr currents were measured in WT and variant CFTR- 
expressing HEK 293T cells, as described in Methods. Panel c. Whole-cell currents of WT-CFTR were measured in HEK 293T cells co-expressed with 
WNK1 and SPAK using patch pipette contained a low concentration of Cr (10 mM). A representative trace of reversal potential measurement is 
shown in the left panel. The permeability ratio Phcos/Pci w3s calculated according to the Goldman-Hodgkin-Katz equation. I-V relationships at the 
indicated points are presented in the accompanying graph. The conductance ratio Ghcos/Gci was calculated by measuring each outward current (i.e., 
slope between Erev and Erev+25 niV). RMP, resting membrane potential. Panel d. Whole-cell currents of R170H-CFTR were measured in HEK 293T cells 
using the same protocol shown in panel c. Panel e. A summary of the Phcos/Pci values obtained from WT-CFTR in the standard state (left) compared 
to WT-CFTR and the nine CF7"R®° variants with WNKl + SPAK activation (right, underlined). Panel f. A summary of the Ghcgs/Gci values in the standard 
state (left) with WNKl + SPAK activation (right). Values throughout are means ± SEM. * p<0.05, **p<0.01: difference from WT in cells co-expressed 
with WNKl and SPAK. 
doi:1 0.1 371/journal.pgen.1 004376.g001 
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Figure 2. Molecular modeling and simulations of CFTR WT and variants. Panels a and b display the side and bottom views, respectively, of 
the WT CFTR dimer, where the two nucleotide-binding domains and the two membrane-spanning domains are labeled as NBD and IVISD. The shaded 
region indicates the location of the lipid bilayer. Color key: black, subunit 1 of CFTR, with residues 1-859; blue, subunit 2, residues 860-1480; red 
CFTR variants studied. Panel c shows the charge distribution around D1152H: this negatively charged residue (left; shown in red space-filling 
representation) is surrounded by several positively charged residues [green), especially on its side of the cavity, creating an attractive force that keeps 
the residue from extending into the cavity. Also shown are other negatively charged residues (red stick or space-filling representation), including 
D385, diametrically opposite to D1 152. Panel d shows the corresponding scene for the variant residue, D1 152H (cyan), which can move toward the 
center of the cavity, thus leading to a constriction in the channel diameter. Channel diameter at the location of variant residues: Panel e 
shows the diameter of the channel at the location of L997, as a function of time, both for the WT (L997, green curve) and the variant (F997, red curve), 
based on closest interatomic distance between L997/F997 and D385. On panel f, the same information for the WT and variant D1 152H is shown. In 
both plots, the pore diameter in the WT is larger than that stabilized in the mutants. The histograms of channel sizes are shown along the right 
ordinate. 

doi:1 0.1 371/journal.pgen.1 004376.g002 



Evidence that CFTR variants are associated with 
pancreatitis 

One of the challenges of genetic association studies is 
determining the eflFect of candidate genetic variants by statistical 
tests when the variant is rare or the mutation effect is uncertain. 
One approach is to increase study power by markedly increasing 
study subject numbers, but this approach is prohibitively expensive 
and not always feasible in rare diseases. Another approach is to 
evaluate the combination of statistical trends linked to studies of 
the functional effects of a variant in a biological system and a 
biologically plausible framework. 

In the current study, 1 1 variants that were previously reported 
to be present in chronic pancreatitis but not CF causing 
[16,17,43-52] underwent functional testing. Only CFTR 
M470V and R1162 (not shown) did not meet criteria of altered 
bicarbonate permeability and/or conductance after WNKl and 
SPAK activation (Figure le-f, discussed below). The remaining 9 



CFTR variants were identified at least twice in pancreatitis 
association studies over the past decade. 

Five variants (R74Q, R75Q, R170H, L967S, and R1162L) 
were located in the hinge region that modulates the collective 
movements of the NBDs with respect to the MSDs (Figure 3). 
R74Q, was previously reported in a single chronic pancreatitis 
patient [53] but not in the CFTR2 database. CFTR R74C) was 
identified by us in two cases and no controls (p = ns) and in one 
case who was a SPINKl N34S carrier (p = 0.006). R75Q, is 
considered to be a non-CF causing mutation according to the 
CFTR2 mutation database [54]. CFTR R75C) was identified in 
61/906 cases and 75/1214 controls (6.3 vs. 6.2%, p = ns) but was 
also detected in nine SPINKl N34S/- mutation carriers (9/55, 
16.4%), with strong combined efiect {SPINKl OR 3.7, SPINK1+ 
R75Q, compound OR 12.2, p 0.002). Two of the nine trans- 
heterozygous cases had been previously reported[18]. R75Q^was 
also identified in four cases with a concurrent severe CF-causing 
mutation and in no compound controls. R170H was first reported 



subunit 1 subunit 2 mode 1 DlOde 2 




Figure 3. Location of selected variant sites with respect to the collective modes of the CFTR. The left panel displays the relative motions 
of individual residues along the collective modes 1 (fop) and 2 (faoffom) intrinsically accessible to the two transporter subunits (residues 71-645 and 
846-1445, respectively). These modes divide the CFTR structure into two groups of residues (colored red and blue in the ribbon diagrams) subject to 
opposite-direction motions. The pale blue/pink indicate the central locations (or hinge regions) mediating the concerted anti-correlated movements 
of the two groups. L967, D1152 and R1162 participate in the hinge region that modulates the concerted anti-correlated (opening/closing) 
movements of the two membrane-spanning domains in mode 2. Residues indicated by filled points (/eft curves) or yellow spheres (right diagrams) are 
R74, R75, R117, R170, L967, L997, D1152, R1162, SI 235 and D1270. 
doi:1 0.1 371/journal.pgen.1 004376.g003 
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in two cases of congenital bilateral aplasia of vas deferens in 
England [53] but is not currendy in the CFTR2 mutation 
database. CFTR R170H was identified in three cases and no 
controls (p = ns). L967S has been reported in a single case of 
azoospermia from the CF mutation database [53] but is not in the 
CFTR2 mutation database. L967S was identified in ten cases (one 
trans-heterozygote), two controls (OR 6.9 p = 0.004), and one 
N34S case carrier. Rl 162L is predicted to be a highly deleterious 
variant by SIFT and damaging by PolyPhen modeling [55] and is 
included in the CFTR2 mutation database and classified as a 
variant not causing CF. Although located in a critical portion of 
the CFTR molecule, the association and functional threshold for 
inclusion as a CFTR'^^' variant were not fully met. 

Two variants (L997F and Dl 152H) appeared to reduce channel 
diameter. L997F is considered a mutation of varying clinical 
consequences for CF, with low rates of pancreatic insufficiency 
and retention of chloride conductance [54] . In this study L997F 
was identified both in the cases (0.7%) and controls (1.0%), 
additionally, L997F was identified in one N34S case carrier and 
three compound heterozygous mutation case carriers, but 
independent statistical association with pancreatitis was not 
demonstrated in this study. D1152H is a mutation of varying 
clinical conserjuence for CF and is associated with low rates of 
pancreatic insufficiency and retention of chloride conductance 
[53]. CFTR Dl 152H was identified in four cases and no controls 
(p = 0.014). Two of these cases were in compound heterozygosity 
with F508del. 

Two variants (SI 235 and D1270N) were on the surface of 
NBD2 (Figure 2). SI 235 is a non-CF-causing mutation[54] and 
was identified in 2.4% of cases and 1.4% of control (p = ns), three 
compound heterozygous cases and one N34S case carrier. While 
this did not reach statistical significance in this cohort, multiple 
previous reports of CFTR S1235R in idiopathic pancreatitis 
patients[56,57] and complex functional features [27] were noted. 
D127()N is of varying clinical consequences for CF, with low rates 
of pancreatic insufficiency and retention of chloride conductance 
[54]. D1270N was identified both in the cases (0.3%) and controls 
(0.2%). Although these variants have been identified in previous 
studies, the effects of these rare variants on altered bicarbonate 
permeability and conductance appear to be weak (Figure 1 e-f) 
and the effect on the function of NBD2 (Figures 2-3) is unclear. 
However, they meet minimal criteria for the class on function 
grounds and contribute to the overall effect on disease risk. 

The final variant (Rl 1711) is located in an extracellular domain 
and has functional effects beyond the other CFTR^^ variants. 
Rl 17H is a complex variant that is associated with CF only when 
found in cis with a T5 tract in intron 8. The CFTR Rl 1 7H variant 
was identified in 22 cases (2.3%) and 8 controls (0.7%) (p = 0.001), 
with only 3 cases and 1 control having the CF-associated 
R117H*T5 haplotype (p = ns), which links the CFTR variant 
R117H to pancreatitis regardless of the intron 8 T5 haplotype. 
R117H*T7/T9 was also identified in 9 of the 80 cases with a 
concurrent severe CF-causing mutation and in no CF carrier 
controls. The R117H variant was the only one with reduced 
chloride current density (Figure lb). While the variant was 
associated with altered bicarbonate permeability and conductance, 
the mechanism is yet to be determined. 

Common CFTR variants previously associated with 
pancreatitis but not confirmed in the current study 

The common polymorphisms M470V, T854T, and Q,1463Cl 
had no significant association with pancreatitis, either individually 
or combined in haplotypes, in contrast to a previous report [58]. 
Haplotypes were determined by counting homozygous carriers of 



each subset (M470V, T854T, P1290P, Q,1463Q, and M470V, 
IVS-T, IVS-TG) and applying Fisher's exact test. The IVS8 T/ 
TG/M470V allele was evaluated in 784 NAPS subjects and 

controls, and no significant associations were found, in contrast to 
another report[59]. The possibility that a series of complex 
haplotypes affect CFTR expression or exon skipping was not 
excluded, but no evidence of direct association was seen in the 
current study or our previous pancreatitis GWAS [24]. 

Thirty-seven of the 81 CFTR variants tested were not 
identified in any cases among the NAPS 2 cohort. The remaining 
variants were also not significantiy overrepresented alone or with 
SPINK 1 or CF mutation carrier. I148T was seen in three cases 
and one control, so an effect could not be detected or excluded; 
the in cis deletion mutation 3199del6 was not detected in any 
I148T carriers. The IVS8T5 variant was identified in 9.9% of 
cases and 8.2% of controls, which is not individually significant. 
There were six N34S/T5 trans-heterozygote controls and no 
cases, but the combined odds ratio (OR 3.9) of the SPINKl 
N34S variant with IVS8T5 was not significantly higher than 
N34S alone. Four additional variants were identified in only one 
patient and no controls: CF mutations 2184delA, 3120+1G>A, 
R1162X and a mutation of varying clinical consequence, 
G1069R. 

Taken together, these genotyping and functional studies provide 
strong rationale for inclusion of nine variants as CFTR^^ class 
members. Although additional variants may be added to the 
CFTR^^ class in the future, the current study did not have the 
very large patient size needed to provide adequate power to detect 
statistically significant changes in additional rare variants. In 
addition, other possible mechanisms of CFTR channel dysfimction 
linked to altered bicarbonate conductance are possible, such as 
mechanisms linked to CFTR R117H. 

Structural significance of CFTR^° mutants 

Structure-based simulations can provide insights into molecular 
driving forces and thereby into the mechanisms of channel 
dysfunction. To better understand the location and structural 
effects of the nine amino acid variants conferring risk of 
pancreatitis and causing dysfunction of the electrophysiological 
response to WNKl-SPAK activation, we developed structural 
models of CFTR and conducted dynamic simulations. Because no 
crystallographic structures for the entire human CFTR are 
currently available, we built a homology model based on the 
structure of a bacterial ABC transporter (Savl866) from Staph- 
ylococcus Aureus [35]. Several computational studies have been 
carried out using models of CFTR and other ABC transporters 
that focus on the structure and/ or gating cycle of the molecule 
and the effect of common mutations/deletions (e.g., F508del in 
CFTR) [4,5,35-37,39,60-64]. Our study is, to our knowledge, 
the first to investigate the multiscale dynamics of CFTR by 
examining both the global motions of the overall protein (with 
ENM) and the local effects of particular CFTR^^ variants (with 
MD). The ENM analysis highlighted the critical positioning of 
R74, R75, R170, L967, and R1162 at the global hinge regions 
(those between the NBD and MSD of transporter in mode 1 , and 
between the two NBDs in mode 2), as evidenced by the 
significant suppression of residue ffuctuations in their close 
neighborhood. Mutations at those sites would thus be expected 
to interfere with the functional dynamics of the channel. Our all- 
atom MD study, on the other hand, showed that a substantial 
constriction could arise in channel diameter with substitutions at 
residues lining the wall of the channel. In particular, the L997F 
and D1152H mutants showed channel pore size reductions in 
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their neighborhoods that would direcdy affect conductance 
properties. 

Defects in CFTR bicarbonate transport 

The fact that all of the pancreatitis-associated variants identified 
by genetic screening in this study resulted in defective WNKl- 
SPAK-activted increase in bicarbonate secretion supports the 
argument that this mechanism is critical for bicarbonate-secreting 
cells that utilizf; CFTR as the primary anion channel. The 
importance of bicarbonate conductance across CFTR at the apical 
membrane is magnified if chloride, but not bicarbonate, conduc- 
tance across the basolateral membrane is minimal, as predicted for 
the pancreatic duct cell [6], since the transcellular anion 
conductance is responsible for fluid secretion. Under basal 
conditions, CFTR-mediated bicarbonate permeability is only 
~20% of chloride, and the capacity for facilitating high 
bicarbonate flux for bicarbonate-secreting tissues is limited. Under 
conditions of low-intraceUular chloride, the WNKl-SPAK path- 
way are activated, and this in turn transforms CFTR into a highly 
bicarbonate-permeable anion channel (Figure 1). The molecular 
mechanisms as to how WNKl-SPAK increases the CFTR 
bicarbonate permeability remain unclear. However, increasing 
evidence suggests that ion permeability of anion channels is not 
fixed and can be dynamically modulated by cellular signaling and 
other events [6,5]. 

The pore of anion channels is believed to ha\-e a large 
polarizable tunnel, where ion selectivity is basically determined by 
the hydration energy of ions and polarizabUity of the channel pore 
[65]. Therefore, in general, the CFTR ion channel is more 
permeable to large anions that are more readily dehydrated [66]. 
However, this cannot be applied to HCO3 . Although the size of 
HCO3 (equivalent radius: 2.1 or 2.43 A) is larger than CI 
(1.81 A), most anion channels, including CFTR, exhibit poor 
HCO3 permeability because of the asymmetrical charge 
distribution of HCO3 [67]. A decrease in the CFTR pore 
diameter, as shown in L997F, can affect the permeability of 
HCO3 in many ways, such as by limiting the accessibility of 
large, asymmetrically charged HCO3 to the channel pore. A 
second mechanism of reducing HCO3 permeability and 
conductance is to inhibit the interaction between CFTR and 
WNKl/SPAK or to reduce the WNKl/SPAK-mediated confor- 
mational change of CFTR. The elucidation of precise molecular 
mechanisms of each mutation will provide insights into the 
understanding of HC03~ conduction in CFTR and also in other 
anion channels. 

Conclusion 

Taken together, these findings support a new class of CFTR 
functional variants with a specific defect in responding to WNKl- 
SPAK activation with increased bicarbonate permeability - 

conductance, dubbed CFTR'^^ . As a class, these nine variants are 
more common in pancreatitis cases than in controls and also had 
evidence of significant risk of pathology in other organs utilizing 
CFTR for bicarbonate secretion. New insight into multiple 
plausible mechanisms were gained by developing a structural 
model of the entire CFTR molecule and by analyzing the 
collective dynamics for wild type and disease-causing variants 
that result in altered channel function. Together, these findings 
provide new understanding of the complexity of pancreatic 
disease related to CFTR-associated duct dysfunction. Identifica- 
tion of members of this new class of CFTR variants on DNA 
sequencing of symptomatic patients in whom a bicarbonate 
channelopathy is suspected may provide insight into disease 



mechanisms and guidance for patient-specific clinical manage- 
ment decisions. 

Materials and Methods 

Study cohort 

The NAPS2 cohort was ascertained, and data were collected as 
described previously [23]. AU patients were prospectively enrolled 
using protocols approved by the appropriate IRBs. Physician- 
confirmed diagnosis of pancreatitis was required for enrollment as 
a case-, wliik; questions on CF, chronic sinusitis, and male infertility 
were included on a case report form administered by a clinical 
research coordinator. DNA and phenotypii: data for patients with 
chronic and recurrent acute pancreatitis (n — 984) and healthy 
unrelated controls (n = 467 from the NAPS2 case-control study 
[23,68] plus DNA from additional healthy controls from 
SomaLogic Inc, (Boulder, CO) (n = 377), the Inflammatory Bowel 
Disease Genetics Consortium (Dr Richard Duerr, University of 
Pittsburgh) (n = 338) and additional University of Pittsburgh 
studies of pancreatitis and pancreatic cancer (Drs David 
Whitcomb and RandaU Brand, University of Pittsburgh) (n = 42) 
[24] were evaluated for a final study cohort of 984 cases and 1224 
unrelated controls. 

Genotyplng 

PRSSl genotyping was done by DNA sequencing [69]. 
SPINKl genotyping was done by sequencing exons 2-3 and the 
flanking regions in a preliminary subset of 745 NAPS2 cases, with 
the entire cohort (cases and controls) genotyped for p.N34S, 
p.P55S and c.27delC using TaqMan assays. The SPINKl C.194+ 
5G>A variant [70] was seen in one patient and one control; 
c.l94+2T>C [71] was not identified in the initial sequencing and 
was not further genotyped. 

CFTR variants for the screening panel were selected from a 
review of published papers and abstracts between 1998 and 2010 
[16,17,43-52] and the open access CFTR mutation database 
based in the Hospital for Sick Children in Toronto (http://www. 
genet.sickkids.on.ca) and Johns Hopkins University (Http:// 
CFTR2.org). 

CFTR genotyping was done using a custom MassARRAY 
iPLEX Gold assay (Sequenom, Inc, San Diego, CA) or custom 
TaqMan Gene Expression Assays (Life Technologies Corporation, 
Carlsbad, CA) through the Genomic and Proteomic Core 
Laboratories at the University of Pittsburgh and verified by 
bidirectional DNA sequencing. All cases and controls were tested 
for each of the 81 selected CFTR variants (Table SI). Variants 
were selected in three stages: the most common CF-causing 
mutations in North America, variations that have been reported in 
pancreatitis literature at least once and a subset of variants that 
have been reported in CF patients but for which the biological and 
pathological relevance remains undetermined (Mutations of 
Undetermined Clinical Significance). 67 SNPs (125GtoC, 
1716G>A, 1717-1G>A, 1898+1G>A, 2183AA>G, 2184delA, 
2789-I-50A, 3120-I-1G>A, 3659delC, 3849-1- lOkbOT, 621-1- 
1G>T, 71 1-I-5G>A, A455E, Dl lOH, Dl 152H, D1270N, D44,SY, 
D579G, F1052V, F1074L, F508C, F508del, G1069R, G1244E, 
G1349D, G178R, G542X, G551D, G551S, I1131L/V, I148T, 
I336K/T, I507del, 1807M, IVS8T5, K1180T, L1065P, L967S, 
L997F, MIV, M470V, M9521, M952T, N13()3K, P67L, 
Q.1463Q, R1070Q, R1162X, R117C, R117H, R170H, R258G, 
R297Q, R31C, R352Q, R553X, R668C, R74W, R75Q, S1235R, 
S1255P, S48,5R, S977F, T338T T854T, V201M, W1282X) were 
multiplexed into 6 wells; 14 SNPs (S492F, S945L, R74Q, R560T, 
R1162L, G85E, I1027T, R334W, R347P, G576A, 711-I-1G>T, 
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1001+1 lOT, P1290P, 3199del6) were ascertained separately via 
TaqMan Gene Expression Assays, with repeat confirmation of all 
positive results. 3199del6 was genotyped via TaqMan on all 
samples that tested positive for I148T. In addition, the intron 8 
boundary was directly sequenced in 873 subjects to determine the 
significance of the IVSB T/TG tract. 

Statistical analysis 

Significant diffc'rc'iicx's in carrier frequencies among cases and 
unrelated controls were determined by chi square analysis or 
Fisher's exact test, and two-tailed p-values are reported. The 
results of each set of experiments are presented as means ± SEM. 
Statistical analysis was performed using Student's ^-tests or analysis 
of variance followed by Tukey's multiple comparison test as 
appropriate. P<0.05 was considered statistically significant. 

Ceil culture and plasmids 

HEK 293T cells were maintained in Dulbecco's modified 
Eagle's medium-HG (Invitrogen, Grand Island, NY) supplement- 
ed with l()"/o fetal bovine serum, 100 U/ml penicillin, and 
0.1 mg/ml str(;ptomycin. The mammalian expressible plasmids 
for hCFTR[72], Myc-rWNKl [21] and Flag-mSPAK [20] were 
described previously. Plasmids were transiendy transfected into 
cells using Lipofectamine 2000 reagents (Invitrogen, Grand Island, 
NY). An average transfection rate over 90% was confirmed by 
transfection with a plasmid expressing green fluorescence protein 
(pEGFP-Nl). Plasmids expressing variant hCFTRs were generated 
using a PCR-based site-directed mutagenesis kit (Stratagene, Santa 
Clara, OA). 

Immunoblotting 

Immunoblotting was performed using conventional methods 
[20]. Briefly, cells were harvested with lysis bufier (20 mM HEPES 
pH 7.4, 150 mM NaCl, 5 mM EDTA, 1% Triton X-100, 1 mM 
NaV04, and 1 mM fi-glycerophosphate) containing a complete 
protease inhibitor mixture (Roche Applied Science, Mannheim, 
Germany). Protein samples were suspended in a sodium dodecyl 
sulfate bulfer and separated by SDS-polyacrylamide gel electro- 
phoresis. The separated proteins were transferred to a nitrocellu- 
lose membrane and blotted with appropriate primary and 
secondary antibodies, and protein bands were detected with 
enhanced chemiluminescence solutions. Antibodies against CFTR 
(M3A7, Millipore, Billerica, MA) and aldolase A (N-15, Santa 
Cruz Biotechnology, Inc., DaUas, TX) were obtained from 
commercial sources. 

Electrophysiology 

Voltage and current clamp experiments were performed on 
HEK 293T cells transiendy transfected with hCFTR as previously 
reported with slight modifications [20]. Briefly, cells were 
transferred into the bath mounted on a stage with an inverted 
microscope (IX-71, Olympus, Osaka, Japan). The pipettes were 
pulled by a Sutter P-57 puller and have free-tip resistances of 
about 2~5 MQ. These were connected to the head stage of a 
patch-clamp amplifier (Axopatch-700B, Molecular Devices, Sun- 
nyvale, CA). Ag-AgCl reference electrodes were connected to the 
bath via a 1.5% agar bridge containing 3 M KCl. Liquid junction 
potentials were corrected for each experimental solution as 
described previously [20]. For the anion permeability test, individ- 
ual data were corrected by measuring the olfset potential shift 
induced by the replacement of anion solution after each 
experiment. The conventional whole-cell clamp was achieved by 
rupturing the patch membrane after forming a gigaseal. Voltage 



and current traces were stored and analyzed using Clampfit v. 10.2 
(Molecular Devices, Sunnyvale, CA). Currents were sampled at 
5 kHz. AH data were low-pass filtered at 1 kHz. 

The high-chloride pipette solution contained (mM) N-methyl D- 
glucamine chloride (NMGD-Cl), 5 ethylene glycol tetraacetic acid, 
1 MgCla, 3 Mg-ATP and 10 4-(2-hydroxyethyl)-l-piperazineetha- 
nesuUbnic acid (HEPES)(pH 7.2). The low-chloride pipette 
solution was prepared by replacing Cl~ with equimolar glutamate. 
The stand bath solution contained (mM) 146 NMDG-Cl, 1 CaCla, 1 
MgClg, 10 glucose and 10 HEPES (pH 7.4). The high-bicarbonate- 
containing bath solution was made by replacing NMDG-Cl with 
equimolar choline-HCOs. The bicarbonate-containing solution was 
continuously gassed with 95% 02-l-5% CO2. 

In all experiments, currents generated by CFTR were 
confirmed by the following three characteristics: 1) activation of 
current by the treatment with cAMP (5 \iM forskolin and 100 |J,M 
l-methyl-3-(2-methylpropyl)-7H-purine-2,6-dione (IBMX), 2) a 
linear I-V relationship and 3) inhibition of current by the 
treatment with the CFTR inhibitor CFTRinh-172. 

The current reversal potential (E^v) was measured either in 
current clamp mode or in voltage clamp experiments. Resting 
membrane potential (RMP) was recorded in zero current clamp 
mode. To test the current-voltage relationship during zero-current 
clamp recording, clamp mode was shifted to the voltage clamp 
mode, and the I-V curve was achieved with ramp pulses from - 
100 to 100 mV (250 ms, holding potential; near tiie RMP). All 
currents were corrected for capacitative currents and the I-V 
relationship was plotted using the values of current density (pA/ 
pF). The relati\'e anion permeability was determined by the 
re\-(-rsal potential shift (AE,-,.,. = Ej„„(X) - Ej,.,.(Cl)) induced by 
replacing extracellular Cl~ with X~ anion using the Goldman- 
Hodgkin-Katz equation as follows: Px/Pci= exp(AErev/(RT/zF) 
- ([Cr]o/[Cr]'„)) _x ([Cr]'y[X-]„), where Pl'^ is the bath 
concentration of CI ; [CI ]„ is the residual CI in the substituted 
solution; [X~]„ is the concentration of substitute ion; and i?, T, z 
and F have their conventional thermodynamic meanings. The 
anion outward chord conductance (Gx: X is anion) between Erev 
and IL^gy +25 mV was achieved by linear plotting. 

Structural modeling 

The homology model of human CFTR (UniProt accession 
code: PI 3569) was obtained using the Swiss-Model Workspace 
software [73]. The most recentiy resolved crystal structure of the 
Staphylococcus Aureus savl866 ABC transporter, fitted to the 
human CFTR (PDB code: 4A82, 2.0 A resolution) [36] was 
adopted as a template, and the structural model was completed 
using the X-ray crystallographic structure of the NMD2 region of 
human CFTR (PDB code: 3GD7, resolution 2.7 A) (see Figure 
S3a). This model de\'iates from the template structure by 1.6 A 
RMSD, and in our simulations the RMSD levels off at ~3.5 A, 
The MSD pore-lining residues and pore radius profile (Figure 
S3b-c) were consistent with those observed in a homology model 
constructed by Norimatsu and coworkers [37,39], which was 
based on an earlier structure (PDB code: 2HYD, resolution 3 A) 
[35]. Using this model for WT CFTR, we generated in silico 
models for tiie mutants L997F and Dl 152H. 

Of note, the collective modes predicted by ANM are highly 
robust and they are not sensitive to small structural variations (like 
those due to a different model). 

IVlolecular dynamics simulations and elastic network 
model analysis 

Molecular dynamics simulations were performed using the 
AMBERl 1 [74] package (GPU version of the pmemd program). 
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with the Amber99SB [75] force field and using the TIP3P water 
model. The protocol consisted of an initial minimization in 
vacuum, using 1,500 steepest descent and 1,500 conjugate 
gradient steps, to remove strong steric contacts, followed by 
another minimization of 5,000 steepest descent and 5,000 
conjugate gradient steps, in explicit solvent, followed by a 
production run of 50 ns. The systems were kept at a temperature 
of 300 K, using Langevin dynamics with a collision frequency of 
2 ps~'; the SHAKE algorithm was adopted to use a 2 fs time step. 
The stability of the system was assessed by verifying the 
convergence of the root mean square deviation (raids) of its heavy 
atoms, after the first 5 ns of simulation. 

As to the pore regions where we examined the local effects of 
substitutions, we allowed for the relaxation and optimization of 
interactions during the described protocol. The simulations, 
thus, gave rise to local rearrangements in the neighborhood of 
the mutation sites and permitted us to extract statistical data on 
the average pore diameter at the constriction zone and its 
fluctuations. 

The elastic network model analysis of collective dynamics was 
performed using the approach reviewed earlier[40]. Collective 
modes of motions are evaluated by eigenvalue decomposition of 
the connectivity/Hessian matrix, using the Gaussian/Anisotropic 
network model. The shape of the mode permits us to identify 
regions subject to large fluctuations as well as domains undergoing 
anti-correlated movements (colored blue and red in the ribbon 
diagrams. Figure 3). 

Ionic diameter 

The radii of the mono-atomic chloride ion was taken from 
Bondi [34]. The equilibrium geometry of bicarbonate ion was 
optimized using ah initio quantum mechanics at DFT level, with 
the B3LYP/6-31 IG** basis set, via the Gaussian 03 software. This 
resulted in a bicarbonate ion that could be fit in a minimum box of 
size 3.40 Ax4.86 Ax5.39 A. This yields a van der Waals radius of 
2.1 (or 2.43) A for the bicarbonate ion, based on the two smaller 
dimensions (or the second largest dimension) that define the 
minimal cross-sec tional area. 

Unless specified otherwise, when we refer to the diameler of the 
pore, we mean the minimal diameter at the specific location of the 
mutation, as opposed to the distribution of diameters along the 
pore. 

Supporting Information 

Figure SI Measurement of Phcos/Pci in cells without WNKl 
and SPAK co-expression. Control experiments were performed. 
Whole-cell recordings were performed to measure CFTR 
bicarbonate permeability by replacing the bath solution with high 
HCO^^-containing (140 mM) solution. Cells were stimulated with 
cAMP (5 nM forskolin and 100 IBMX) after establishing 
whole-cell configuration. The current to voltage relationship (I/V 
curve) was obtained by depolarizing ramp pulses from — 100 to + 
100 mV. The permeability ratio Phcos/Pci was calculated 
according to the Goldman-Hodgkin-Katz equation. I-V relation- 
ships at the indicated points are presented in the right panel. The 
conductance ratio Ghco^/Gci was calculated by measuring each 
outward current (slope between E„.v and Ej„.+25 mV). Replacing 
the bath solution with a high-bicarbonate-containing solution 
induced a strong positive shift in Erev, indicating that bicarbonate 
permeability is much smaller than that of chloride. Summarized 
results of multiple experiments (n = 7) are presented in Figure le 
and f. 
(PDF) 



Figure S2 Representative current to voltage (I-V) plots otCFTR 
variants in whole-ceU current measurements. Whole-ceU record- 
ings were performed to measure the CFTR HCO3' permeability 

and conductance by replacing the bath solution with high HCO3"- 
containing (140 mM) solution. The pipette solution contained 
10 mM Cr. WNKl and SPAK kinases were coexpressed with 
wild-type (WT) or variant CFTR. Cells were stimulated with 
cAMP (5 ixM forskolin and 100 jxM IBMX) after establishing 
whole-cell configuration. The I-V curve was obtained by 
depolarizing ramp pulses from —100 to +100 mV (250 ms), and 
all [:urrents were txjrrected for capacitative currents. Treatment 
with the CFTR inhibitor CFTRi„h-172 (20 nM) inhibited the 
HCO3" currents by an average of 91.8±3.0% (WT-CFTR with 
WNKl & SPAK coexpression, measured at +100 mV, n = 4) 
indicating that CFTR mediates most of the HCO3' currents. The 
permeability ratio Phc:03/Pci was calculated according to the 
Goldman-Hodgkin-Katz equation. The conductance ratio 
Ghcos/Gci was calculated by measuring each outward current 
(slope between Erev and Ejev-i-25 mV). 
(PDF) 

Figure S3 Structural features of the present homology model 
and comparison with previous work. Panel a shows the 
superposition of the current homology model onto the X-ray 

structure of ABC transporter from S. aureus, which yields an 
RMSD of 1.6 A. Panel b shows the position of some residues 
lining the pore at the MSD, located on two TM helices, 6 (red) and 
12 (yellow). These residues were reported by Sansom and 
collaborators to be lining the MSD pore [61], based on both, a 
homology model and experimental cystidnc scanning. Panel c 
shows the MSD pore radius profile — measurc'd as the radius of the 
smallest circle that fits the cross sectional area at each elevation 
along the pore axis (perpendicular to the membrane plane), and its 
extension toward the cytoplasmic region. The MSD portion 
(indicated by the upper abscissa label) is comparable to that 
reported earlier by Norimatsu et al [37]. 
(TIFF) 

Figure S4 Contribution of the slowest modes to the square 
displacements of residues in CFTR. Square displacements are 
calculated using the ENM representation of the two subunits of the 
structure (residues 71-645 and 846-1445, respectively). The right 
panel shows a color-coded ribbon diagram where regions subject 
to large fluctuations are colored pink, and those maximally 
constrained, blue. Note that R74, R75, R170, L967, D1152 and 
Rl 162 lie in the highly constrained region. 
(TIF) 

Table SI 81 CFTR variants genotyped in pancreatitis patients. 
The CFTR mutations investigated in this study are reported with 
legacy nomenclature and relative ranking among the American 
College of Medical Genetics most common classic cystic fibrosis- 
causing mutations found in North America (CF). Those found to 
be associated with cases in the current cohort include an X in the 
Pane Disease column. *IVS8 T5 and Rl 17H are reported but CF 
disease causing only when in cis with each other or 1VS8 T5 with 
IVS8 TG12orl3. Intronic mutations are reported in standard 
nomenclature "####+/-##N>N" except IVS8-T5, (1210- 
12T[5]). 
(DOCX) 
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